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Abstract 

We review the fundamentals of coupling constant metamorphosis 
(CCM) and the Stackel transform, and apply them to map integrable 
and superintegrable systems of all orders into other such systems on 
different manifolds. In general, CCM does not preserve the order of 
constants of the motion or even take polynomials in the momenta to 
polynomials in the momenta. We study specializations of these actions 
which do preserve polynomials and also the structure of the symmetry 
algebras in both the classical and quantum cases. We give several 
examples of non-constant curvature 3rd and 4th order superintegrable 
systems in 2 space dimensions obtained via CCM, with some details 
on the structure of the symmetry algebras preserved by the transform 
action. 

PACS: 02.00.00, 02.20. Qs,02.30.Ik, 03.65.Fd 
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1 Introduction 



There has been a recent rapid expansion in the number of known classical 
and quantum superintegrable systems of order 2, [T],[2], and, particularly, of 
order 3 and higher, O SI [5], El [3 E] • For many of these systems it has been 
demonstrated that the algebra generated by the fundamental higher order 
symmetries closes under the Poisson bracket in the classical case, and under 
the commutator in the quantum case, to form a finite dimensional quadratic 
or cubic algebra. The representation theory of these algebras and their asso- 
ciation with basic properties of the special functions of mathematical physics 
is of great current interest, [U [TUl [TT], [T21 [ISl El [IS]- Indeed, the basic prop- 
erties of Gaussian hypergeometric functions and their various limiting cases, 
as well as Lame, Mathieu and Heun functions, and ellipsoidal harmonics all 
appear as associated with 2nd order superintegrable quantum systems via 
separation of variables. These functions as well as orthogonal polynomials of 
a discrete variable, including the general Wilson and Racah polynomials, are 
bound up with function space models of the irreducible representations of the 
quadratic algebras associated with 2nd order superintegrable quantum sys- 
tems. The Painleve transcendents (not associated with variable separability) 
appear in the study of 3rd order superintegrable systems. Some examples 
are known for conformally fiat manifolds in n dimensions, piUl [T71 [TH] , but 
most results are known for 2 and 3 dimensional conformally flat spaces. 

There is a disconnect, however, between what is known for 2nd order 
superintegrable systems and what is known for 3rd and higher order sys- 
tems. For 2nd order superintegrable systems, classical and quantum, all 
such systems and all manifolds of dimension 2 on which they occur have 
been classified and the mechanism of the closure of the quadratic algebra 
is well understood, [201 [211 [221 [23] • For conformally fiat manifolds in 3 
dimensions great advances have been made although the classification and 
structure analysis is not yet complete, pH |25l [211 [261 [27]. A major tool for 
obtaining these 2nd order results has been the Stackel transform [28] , a vari- 
ant of coupling constant metamorphosis [22], which enables a 1-1 invertible 
transformation between a 2nd order superintegrable system on one manifold 
and a superintegrable system on another manifold that preserves the symme- 
try algebra structure. This has given us an elegant method for classification 
of all 2D superintegrable systems through the important fact that every such 
system can be shown to be the Stackel transform of a system on a constant 
curvature space, [21 [201 [23] ■ Also it gives important insight into the structure 
of Koenigs' remarkable potential-free results [30] . Similar results are known 
for 3D systems but the classification is not yet complete, [251 126l [27] . 

For 3rd and higher order superintegrable systems, however, there is no 
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structure and classification theory. Only examples are known, and these are 
very difficult to obtain. The symmetry algebras can be computed for each 
example but the mechanism for their closure and structure is not understood. 
Virtually all known examples are in 2D or 3D Euclidean space. The present 
paper is a first attempt at refining a tool (CCM/Stackel transform), that 
has proved so successful in the classification and structure theory for 2nd 
order systems, so that it applies to higher order superintegrable systems. 
There are two basic issues here. The first is that COM in general doesn't 
preserve the structure of the symmetry algebras. We have to determine a 
suitable restriction that does preserve the structure. Secondly, CCM is a 
classical phenomenon; its extension to the quantum case is not automatic 
and requires special care. In this paper most of our classical results will be 
stated for n dimensional systems whereas, for simplicity, the quantum results 
will be limited mostly to 2 dimensions. 

In future papers we will extend the operator CCM to 3 and higher dimen- 
sions and employ this tool to attack the structure and classification theory 
for 3rd and higher order superintegrable systems in all dimensions. An im- 
mediate result of the present paper is the explicit display of a large number of 
higher order superintegrable systems on manifolds not of constant curvature, 
the existence of which seems not to be widely recognized. We also provide 
new examples of explicit structure computations for the quadratic algebras 
of some 3rd and 4th order superintegrable systems on 2D Euchdean space 
that map to isomorphic systems on nonconstant curvature spaces. 

Before proceeding to our results we give some basic definitions that we 
employ throughout the paper. A classical superintegrable system on an n- 
dimensional real or complex Riemannian or pseudo-Riemannian manifold is 
defined by its associated Hamiltonian function H = Yliij Q^^PiPj + ^(^) 
the phase space of this manifold. Here g'^^{yi) is the contravariant metric 
tensor in local coordinates x and F(x) is a prescribed function that may de- 
pend on some parameters. The system is superintegrable if it admits 2n— 1 
functionally independent generalized symmetries (or constants of the motion) 
iSfe, /c = 1, • • • , 2n — 1 with Si = 7i where the Sk are polynomials in the mo- 
menta p^-. That is, {n,Sk} = where {f,g} = YTj=i{dxJdp^g - dpjd^^g) is 
the Poisson bracket for functions /(x, p), g'(x, p) on phase space. It is easy to 
see that 2n — 1 is the maximum possible number of functionally independent 
symmetries and, locally, such (in general nonpolynomial) symmetries always 
exist. Most authors, but not us, also demand that the system is integrable, 
i.e., there is a subset of n functionally independent polynomial symmetries, 
say (Si, • • • , Sni such that {Sj, S(} = 0, 1 < s, ^ < n. If the maximum order 
of the polynomials corresponding to the generating symmetries is N , we say 
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that the system is iVth order superintegrable. 

Superintegrable systems can lay claim to be the most symmetric Hamilto- 
nian systems though many such systems admit no group symmetry; the sym- 
metry is "hidden". Gcncrically, every geometrical trajectory in phase space 
(but not the time dependence of the trajectory p(t),x(t)) of the Hamilton 
equations of motion for the system, is obtained as the common intersec- 
tion of the (constants of the motion) hypersurfaces »Sfc(p, x) — Ck, k — 
0, • • • , 2n — 2. The orbits can be found without solving the equations of 
motion. Since every known superintegrable system is also integrable, this is 
better than integrability. A case can be made that the 2nd order supcrinte- 
grability of the Kepler- Coulomb two-body problem, forcing the existence of 
conic sections as trajectories, is the reason that Kepler was able to determine 
the planetary elliptical orbits before the invention of calculus. 

There is an analogous definition of superintegrability for quantum systems 
with Schrodinger operator 



the Laplace-Beltrami operator plus a potential function. Here it is required 
that there are 2n — 1 functionally independent differential operators, 5*1 = 
if, 5*2, • • • , S2n-i such that and [H, Sk] = HSk — S^H — 0. Often there is 
a 1-1 relationship between classical and quantum superintegrable systems 
associated with a potential and then functional independence refers to the 
classical system. In those cases where there is no classical analog, however, 
there is no agreed upon definition of quantum functional independence. A 
basic motivation for studying these systems is that they can be solved ex- 
plicitly, often in multiple ways. Typically their symmetry algebras close to 
form quadratic, cubic, or similar algebras whose representation theory yields 
spectral information about the quantum system. 

In the following sections we review the basic definition and properties 
of couphng constant metamorphosis (CCM) and the closely related Stackel 
transform for classical systems. These concepts apply to any Hamiltonian 
system with potential, not just superintegrable systems. Then we define 
specializations of these general concepts that preserve the order of symmetries 
and also define symmetry algebra isomorphisms. It is these specializations 
that are needed for the study of superintegrable systems. Then, and most 
importantly, we find quantum analogs of these classical transforms. At each 
stage we provide examples, several of them new. 
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2 Coupling constant metamorphosis 



The basic tool that we will employ follows from "coupling constant meta- 
morphosis" (CCM), a general fact about Hamiltonian systems, pointed out 
in [29]. Let 7i(x, p) + aU{x.) define a Hamiltonian system in 2n dimensional 
phase space, with canonical coordinates Xj,pj. Thus the Hamilton- Jacobi 
equation would take the form 7i(x, p) + af/(x) = E. Assume that for every 
value of the parameter a the system admits a constant of the motion /C(a), 
analytic in a. 

Theorem 1 Coupling constant metamorphosis. The Hamiltonian Ti' = (7Y— 
E)/U admits the constant of the motion IC' = /C(— 7^'), where now E is a 
parameter. 

PROOF: Note that if F, G are functions on phase space of the form G(x, p), 
F = F[a) = F{a, x, p) where a = q;(x, p) then 

{F,G} = {F(a),G}U=«(x,p) + 9,F(a)U=„(x,p){a,G}. 

By assumption, {/C(a),7i} = — a{/C(a),f/} for any value of the parameter 
a. Thus 



Now 



{K:{n'),n'} 



dMa){n', n'} + iEiJ^^(n' + «) 



= 0. 

a=-H' 



Q.E.D 



Corollary 1 Let ICi{a) , IC2{a) be constants of the motion for the system 
7i(x, p) +a;f/(x). Then {/Ci,/C2}(tt) = {A^i(tt), A^2(a)} is also a constant of 
the motion and 

{ic,i~n'),}C2i-n')} = {/Ci,/C2}(-7^')- 

Clearly CCM takes integrable systems to integrable systems and super- 
integrable systems to superintegrable systems. We are concerned with the 
case where 

n 

n=J2 + ^(^) + "f^(^) ='Ho + V + aU (1) 
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is a classical Hamiltonian system on an n-dimensional pseudo-Riemannian 
manifold and are interested only in those constants of the motion /C that 
are polynomial in the momenta. As we shall see, in the case of 2nd order 
constants of the motion there is special structure. The 2nd order constants 
of the motion are typically at most linear in a, so they transform to 2nd 
order symmetries again. In this case CCM agrees with the Stackel transform 
that we shall take up in the next section. However, in general the order of 
constants of the motion is not preserved by coupling constant metamorphosis. 

Example 1 The system 

H = pI+pI + bl^/x^+ b2X2 

admits the 2nd order constant of the motion /C^^^ = + ^23^2 ond the 3rd 
order constant of the motion IC^^^ = pl + ^biy/xipi — ^p2, (14^ '^'^^ references 
contained therein). If we choose all = a^/xi then the transform of IC^^^ 
will be 5th order. If we choose all = ax2 then the transform of IC^^'' will be 
rational, but nonpolynomial. Thus to obtain useful structure results from this 
general transform, and to obtain results that have the possibility of carrying 
over to the quantum case, we need to restrict the generality of the transform 
action. 

3 The Jacobi transform 

Here we study a specialization of coupling constant metamorphosis to the 
case where = 0. The special version of the transform we study takes 
iVth order constants of the motion for Hamiltonian systems to A^th order 
constants of the motion. An A^th order constant of the motion /C(x, p) for 
the system 

n 

^ = I] 9''P^PJ + f/(x) =no + U (2) 

is a function on the phase space such that {JC,H} = where 
/C = /Cat + }Cn^2 + }^N-4: + • ■ ■ + JCq, n even, 

/C = /Cat + 1Cn^2 + + 1- ^1, n odd. 

Here, /Cat 7^ and ICj is homogeneous in p of order j. This implies the 
conditions 

{/C^,7^o} = 0, (3) 
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{/C^_2fc, U} + {/C^_2fc-2, no} = 0, k = 0,l,--- , [N/2] - 1, (4) 
and, for odd, 

{/Ci,f/} = 0. (5) 
The case = 1 is very special. Then /C = /Ci and the conditions are 

{/C,7^o} = 0, {/C,t/} = 0, 

so /C is a Kilhng vector and U is invariant under the local group action 
generated by the Killing vector. 

For N = 2, K. = }C2 + JCq and the conditions are 

{/C2,7^o} = o, {/C2,t/} + {/Co,7^o} = o, (6) 

so /C2 is a 2nd order Killing tensor and U satisfies (linear) Bertrand-Darboux 
integrability conditions. 

For N = 3, K. = K.3 + K-i and the conditions are 

{/C3, Ho} = 0, {/C3, U} + {}C^,no} = 0, {/Ci, f/} = 

so /C3 is a 3rd order Killing tensor. The integrability conditions for the last 
2 equations lead to nonlinear PDEs for U. 

Theorem 2 Suppose the system ^ admits an Nth order constant of the 
motion K, where N >1. Then 

^ = J2( — 77 — ) ^^-2i 

is an Nth order constant of the motion for the system (Tig — E)/U. 
PROOF: It follows from the general conditions ®, (HD, ©• That 

[7V/2] 

/C(a) = ^ a^lCN~2j 

j=0 

is a constant of the motion for the system T-C{a) = TYq + c^U. Then from 
Theorem [H we have that K,{— ^°~^ ) is an iVth order constant of the motion 
for the system {Hq - E)/U. Q.E.D. ^ 

Note that if we set = then K, becomes an A^th order Killing tensor 
for the free system TCq/U. 
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Corollary 2 Suppose the system Ti.o + U is Nth order superintegrable. Then 
the free system TCq/U is also Nth order superintegrable. 

We will call /C a Jacobi transform of /C, in recognition of its close rela- 
tionship to the Jacobi metric, [31] page 172, and to distinguish it from the 
Stackel transform and more general coupling constant metamorphosis. Note 
that the Jacobi transform for general parameter E is invertible. 

Corollary [2] tells us that each of the 3rd order superintegrable systems 
found by Gravel in 2D Euclidean space, [4], yields superintegrable systems 
on conformally flat manifolds, usually not of constant curvature. 

Corollary 3 The Jacobi transform satisfies the properties 

{^} = {IC,C}, KX = tt, 

and, if JC, C are of the same order, alC + bC = aK, + bC. Thus it defines a 
homomorphism from the graded symmetry algebra of the system TYq + U to 
the graded symmetry algebra of the system {Hq — E)/U . 

Example 2 Consider the system of Example\^ = Pi+pi + bi^/xi + b2X2, 
and let U = biy^ + b2X2 + b-s for some fixed bi, 62, &3 with 6162 7^ 0. The new 
Hamiltonian is 

^ pI + pI-e 

^=b,^, + b2X2 + b, 

and the Jacobi transforms of }C^'^\ JC^^^ are 

^ ^^^^ (&iv/^+ ^23^2 + ^3) ' 

2 462 V^iv^ + "2a;2 + 03 / 



4 The Stackel transform 

Using the same notation as in the previous section, and a particular nonzero 
potential U = V{x, bo) we define the Stackel transform for a system TC = 
T^o + l^(x, b), p8]. The transform of /C = /Ci is /C^ = /Ci The transform of 

/C = /C2 + /Co is /C = /C if-T^- (Here ]Cj is a homogeneous polynomial 

in p of order 2j, and K.J is the restriction of Kj to the potential V = U.) 
The transform maps 1st and 2nd order constants of the motion for Ti to 
constants of the motion for the system Ti/U. Thus the system Ti is 2nd 
order superintegrable iff the system Ti/U is 2nd order superintegrable. For 
completeness we review briefly the direct proofs of the basic theoretic facts. 
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Theorem 3 Let K, he a 2nd order constant of the motion for the system Ti 
and U he a particular instance of the potential V . Then IC is a 2nd order 
constant of the motion for the system Ti/U . 



PROOF: 



^} = - ^} = ({^2, U] + {/Co^ Ho}) = 



Q.E.D. 



Corollary 4 Let IC, C he 2nd order constants of the motion for the system 
H and let C he their respective Stdckel transforms determined hy the 
potential U. If {/C, C} = then {JC, C} = 0. 



PROOF: Suppose {/C, C} = 0. We have 

{/C, £} = {/C2, £2} + ({/C2, Co} + {/Co, C2}) = 0, 

where the first term on the right is of order 3 and the second term is of order 
1. Thus 

{/C2,£2} = {/C2,/:o} + {/Co,/:2} = o. 

Then, a straightforward computation yields 

{IC, C} = {IC, £} - ^ ({/C2, £0^} + {/Co^, £2}) = 0. 

Q.E.D. 

Corollary 5 Let {IC, C} = he as in Corollary ^ and assume that one 
instance of the potential V is the constant 1, i.e., ^(bi) = 1. Then if 
{IC, C} = we must have {IC, C} = 0. 



PROOF: Suppose {IC, C} = 0. Then the order 3 and order 1 terms on the 
left hand side of this expression must vanish separately: 

{IC,, £2} - ^ ({/C2, ^o''} + {/Co^ A}) = 0, (7) 

{IC2, Co} + {/Co, C^} - ^ ({/C2, + {/Co'', C2}) = 0, (8) 
Identity (JTj) says that 

{iC2,C2} = noX (9) 
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where X = ^ (^{)C2, Cq} + {JCq , C2}) ■ Since C2 are 2nd order Killing 
tensors of it follows easily from the Jacobi relation for the Poisson bracket 
that A" is a Killing vector. From identity (E]) we obtain the result 

{1C,C} + 'HX = {iCX} = (10) 

Taking the Poisson bracket of the left hand side of this last identity with Ti 
we see that A* is a first order constant of the motion for system Ti. From ([8]) 
we have 

for any nonzero choice of potential V . Choosing \^ = 1 we find 

X = {lC2,Cl} + {lClC2}. (11) 

From relation (ED with V" = 1 we have {/C2, 1} + {/Cj, Hq] = so {/C^, Hq} = 
0. Since the metric is nondegenerate, this implies that /Cg = ci, a constant. 
Similarly, £q = C2 is constant. Thus (|TT|) implies X = 0, which together with 
(dUD implies {IC, C} = 0. Q.E.D. 

An alternate way of proving Corollary [5] is to demonstrate that there is an 
"inverse" Stackel transform that takes the system Ti/U to Ti via the special 
potential 1/U . The outcome of applying the initial transform to a 2nd order 
constant of the motion /C of 7i and then transforming back is /C — /Cq^, 
where /Cq is a constant. If each 2nd order symmetry /C is normalized by the 
requirement /Cq = 0, (by adding a suitable constant) then this action is the 
identity operator. 

These results show that the Stackel transform takes 2nd order superin- 
tegrable systems to 2nd order superintegrable systems, preserves variable 
separability, and is invertible. As stated in this generality for second order 
symmetries, the Stackel transform is not a special case of coupling constant 
metamorphosis, although the two transforms are closely related. However in 
the situation where the potential functions l^(x, b) form a finite dimensional 
vector space, which is usual in the study of 2nd order superintegrability, then 
the transforms coincide. In this case, by redefining parameters if necessary, 
we can assume V is linear in b. 

Now we will investigate extensions of the Stackel transform to higher 
order constants of the motion, under the assumption that V^(x, b) is linear 
in b = (60, ^m), U is of the form f/(x) = V^(x, b°) and the potentials 
V^(x, b) span a space of dimension M + 1. In particular, 

M 

l^(x,b) = 6o + l]f/^*nx)&. (12) 

i=l 
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where the set of functions {1, U^^\x), ■ ■ ■ , U^^^^x.)} is hnearly independent. 
In the study of 2nd order superintegrabihty, typically the 2nd order constants 
of the motion are linear in the b and the algebra generated by these sym- 
metries via products and commutators has the property that a constant of 
the motion of order N depends polynomially on the parameters with order 
< [iV/2]. Thus we consider only those higher order constants of the motion 
of order N of the form 

[N/2] 

/C = ^ /CAr_2i(p, b) (13) 
j=0 

where /CAr„2i(a'P, b) = a^^'^^ ICN-2j{p,^) and /CAr_2j(p, ab) = a^)CN-2j{p,^) 
for any parameter a. Let /C(b) be such an A^th order constant of the motion. 
Then 

/C(a) = /C(p,b + ab(°)) (14) 

is an A^th order constant of the motion for the system with Hamiltonian 
T^o + ^(x, b) + aUix). Applying Theorem [1] we have 

Theorem 4 Let /C be an Nth order constant of the motion for the system 
Ho + V{yi, b) where V is of the form ( fl^) and /C is of the form ( fl^) . Let 
]C{a) be defined by ( [7^ . Then 

is an Nth order constant of the motion for the system (Tio + ^(x, b))/f/(x), 
where 

iCN-2j{ap,h) = a^~'^^ICN-2jiP,h), ^Af_2j(p,ab) = }CN-2j{p,h) (15) 
for any parameter a. 

Example 3 This example of a 4th order superintegrable system is taken 
from 18] and corresponds to the choice k = 2 for the potential V = Ar"^ + 
B/r"^ cos^ (kt) + C/r'^ sin^(A;t) for suitable A, B, C , as written in polar coordi- 
nates. The structure relations and transform are new. Let 

^ = Pi + P2 + + ^2) + ^JZa—ayi + ^ ^2^2 
There are two basic constants of the motion, one of 2nd order, 

22 / 4 I 4\ 

Olj Oil ( 3y ~\ Ou ) 

/C2 = {xm - xmf + y_l + c'-^-^ 
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and one of 4ih order, 



-2c 





{IC2, IC4} .The relations are 

{/C2,7^} = 32(7^2_2/C4)/C2-64(6+2c)/C4+64(6-c)7^2_l28a^)K:l-128a6(6+2c), 
{/C4,7^} = 32/C4(/C4-:^^) + 128a/C27^^-384a2/C^ + 128a6/C4-64(6 + 4c)a7^^ 



7^2 = 64IC2lCiin^-ICi)-64bn^+128{b~c)ICin^-64{b+2c)}Cl-128alCln^+ 

256a2/C^ - 256a6/C2/C4 + 128a(6 + 4c)7^2/C2 + 2560^(6 - c)/C^ - 256a6(6 + 2c)/C4 
+256a(76c+6^-2c^)7^2_256a2(62+4c2+206c)/C2-256a2(2c+6)(62+166c-4c2). 



is a/so superintegrable with 4th and 2nd order generating constants of the 
motion. Note that for B = C = 0, D = 4, the transformed system is defined 
on a Darboux space of type 3, whereas ifB = D = 0,C = l the transformed 
system is defined on a Darboux space of type 2, J^. 

5 2D quantum symmetries 

Here we begin the study of quantum symmetries. The quantization is much 
simpler in the 2D case than for dimensions greater than 2, and for 1st and 2nd 
order symmetries, so we begin with these special cases to gain insight. Here 
the metric, expressed in Cartesian-like coordinates, is ds"^ = X{x){dxi+dx2), 
and the Hamiltonian system Ti. = {pi +P2)/X{^) + V{x.) is replaced by the 
Hamiltonian (Schrodinger) operator with potential 



+256a^(2c - 6)/C2 + 12802(6^ + 40c^ + 206c) 



There is a Casimir constraint 




+ d 



1 



{du+d22) + V{^). 



(16) 



H 



A(x) 
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5.1 2nd order operator symmetries 

A 2nd order symmetry of the Hamiltonian system AT = Ylt j=i {^)PkPj + 
H^(x), with a''^ = a^^ , corresponds to the operator 

1 2 

= V a, (A(x)a^'^(x)5,) + iy(x), a^^ = a^K 
^^"^^ k,j=i 

These operators are formally self-adjoint with respect to the bilinear product 
<f,9>\^ J /(x)c/(x)A(x)dxidx2 

on the manifold, i.e., 

< /, Hg >A=< Hf, g >x, < f, Kg >a=< Kf, g >x 

for all local C°° functions /, g with compact support on the manifold, where 
the domain of integration is or TZ^. If the functions defining a differential 
operator arc singular on a 1-dimensional or 0-dimensional set, we we restrict 
the support of f,g to be bounded away from this set. We define the formal 
adjoint T* of a linear operator T on the space by 

<T*f,g>,=<f,Tg>, (17) 

for all f,g ^ C^. The operators H,K are formally self-adjoint: H* — 

H, K* = K. 

If the Schrodinger equation admits a multiplicative separable solution in 
particular coordinates xi,X2 then the Schrodinger operator can be written 
as 

" = XC)(x.)!x.^)(..) + + ^'"(^') + ^'''(^^)) ^''^ 
where the 2nd order symmetry responsible for the separation is 

+ X(2)(x2)\/«(a;i) - X^'\x,)V^^\x2)) . 
Thus the metric is A(x) = X^^\xi) -\- X^'^\x2) and the potential is V^(x) = 

A 1st order symmetry of the Hamiltonian system C — Y^k=i^^^'^)Pk 
corresponds to the operator 

' / a,(A(x)a'=(x)) 



fe=l ^ 



2A(x) 
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It is easy to show that L is formally skew-adjoint, i.e., L* = —L. 

The following results that relate the operator commutator [A, B] — AB — 
BA and the Poisson bracket are straightforward to verify. 

Lemma 1 

{n, /C} = <^ [H, K] = 0. 
This result is not generally true for higher dimensional manifolds. 
Lemma 2 

{n,C} = ^ [H,L]^0. 

The classical Stackel transform for these systems can easily be extended to 

the operator case. Suppose \^ is a parametrized potential and let U he a, 
special instance of that potential. Let K = di{Xa^^dj) + W = K2 + Kq, 
where Kq = W, be a 2nd order formally self-adjoint symmetry operator of 
H and Kq be the restriction of Kq to V = U. Then 

is the corresponding formally self-adjoint symmetry operator oi H — U~^H, 
with respect to the metric ds"^ = U\{dxl+dxl). Here the order of operators in 
a product is important and a function represents the operation of multiplying 
on the left by that function. 

Theorem 5 1. 

[H, k]^0 ^ [H, K] = 0. 

PROOF: 

1. This is a straightforward verification, using the identities 

[i/o,i^2] = 0, [H^ + V,K2 + K^]=Q, [Ho + U,K2 + K^]=Q, (20) 

[A BC\ = B[A, C] + [A, B]C, [A, C/"^] = -U'^IA, Up'^ 
for linear operators A,B,C and nonzero function U. 

2. This follows from the fact that diKl^ = diW^ = A a'^djU. 
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Q.E.D. 

Note that the second part of the theorem shows that K is indeed for- 
mally self-adjoint on the manifold with metric U\{dx\ + dx^). Another way 
to see this is to use the formal definition of adjoint. With respect to the 
inner product on the space with weight function UX we have < f,g >u\=< 
f,Ug >A=< Uf,g >A. Thus 

< kf, g >ux=< {UK - UK^U~'H)f, g >a=< /, {KU - HK^)g >a 

=< /, {U~^KU - U-'HK^)g >ux • 

This shows that K* = U'^KU - U'^HK^ = K, where the final equality 
follows directly from identities ( !20|) . 

Corollary 6 If K^^\ K^"^^ are 2nd order symmetry operators for H , then 

Since one can always add a constant to a potential, it follows that 1/U 
defines an inverse Stackel transform of H to H. 



5.2 3rd order operator symmetries 

A classical 3rd order symmetry takes the form /C = /C3 + /Ci where 

2 2 

k,j,i=l 1=1 

The conditions {K^.Tio] = are 

2af = -3 ((lnA),a*** + {\nX)ja^'') , z ^ j (21) 
3af' + af = -3 ((In \)ia''' + (In X)ja'^^) , t^j 
2(aP + an ^ -(lnA)iai2^-(lnA)iai"-(lnA)2a'''-(lnA)2ai^2^ 

which are just the requirements that the a'^-^' be the components of a 3rd 
order Killing tensor. The conditions {/C3, V} + {/Ci,7io} = are 

2 

6^ + 52 = 3 ^ Aa^' V„ (22) 

s=l 

2 2 



2 ^ 2 

s=l s=l 
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and the condition {/Ci, = is 



2 



(23) 



s=l 



Now let's consider a 3rd order operator symmetry K that is skew-adjoint. 
The detailed conditions [K, H] = are complicated. However, we will re- 
strict ourselves to systems with potentials that simultaneously admit a 3rd 
order classical symmetry and the related 3rd order quantum symmetry. A 
characteristic feature of such systems, and one that we will exploit, is that if 
U is such a potential then so is all for all scalars a. If /C is the classical sym- 
metry then we can write the operator symmetry in the form K = K'^ + K[ 
where K^,K[ are skew-adjoint of respective orders 3 and 1, 



and the terms a'^-'* satisfy fl^T]) . 

Now replace V by aU. Then the symmetry condition is [K{a), Ho+aU] = 
for all a where K = K'^ + K[{a). We assume that K{a) is analytic in a 
about a = 0. Then, K'^ is independent of a and the dependence of K[ on a is 
at most 1st order. Thus we can write B^{a) = d + aU or K[{a) = K'( + aKi. 
The symmetry condition can be written as 



Note that the 2nd order terms in the 2nd identity are precisely the classical 
conditions f l22l) . The 3rd identity is precisely the classical condition ( l23ll . 
The operator K'( determines the transition from the classical constant of the 
motion to the operator symmetry. Now define 





= + K'l + aKi, Ho + aU] 
= [K', + K'l, Ho] + a{[K', + K'l, U] + [K^, Ho]) + a^[Ki, U] 



for all a. Setting = + K'^ we have the identities 

[Ks, Ho] = 0, [Ks, U] + [K,, Ho] = 0, [K,, U] = 0. 



(24) 



K = K,-K,U-\Ho + b), 



where the operator order is important and 6 is a constant.. A straightforward 
computation using identities, yields [K, U~^{Ho + b)] = 0, so is a 3rd 
symmetry operator for the Hamiltonian U~^{Ho + b). 
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Theorem 6 Let H{a) = Hq + all, let K{a) be a 3rd order skew- adjoint 
symmetry of H , analytic in a about a = 0. Then there are 1st and 3rd order 
skew-adjoint operators Ki,K^ such that K{a) = + aKi and identities 
l[24\ ) are satisfied. The operator K = K3 - K\^U-\Ho + b) is a 3rd order 
symmetry for the system H = U~^{Hq + b). 

Corollary 7 K* = —K so K is a 3rd order formally skew-adjoint symmetry 
ofH. 

PROOF: This is a consequence of K* = —K,H* = H and relations fl^ . 
Q.E.D. 

Note: The preceding argument has to be modified in the special case that 
the system admits a 1st order a-independent symmetry L: [L, H{a)] = 0. 
Then K[{a) need not be at most 1st order as a polynomial in a. Indeed we 
can add a term f{a)L to K[ without changing the commutation relations. 
However, the conclusion fl24p remains correct. 

Example 4 (The 9-1 anisotropic oscillator) Let H{a) = du -\-d22 + a(9x^ + 
Xg). This is a superintegrable system with generating 2nd and 3rd order 
symmetries 

1 = 822 + axl, K = {xid2 - X2di, 822} + -{{xl, di} - 9{xixl, 82}), 

where {Si, S2} = S1S2 + S2S1. Let U = {dxf + x^) + c. It follows that the 
system 

is superintegrable with one 2nd and one 3rd order symmetry. 

5.3 4th order operator symmetries 

Next we consider the case of a 4th order constant of the motion 
2 2 
}C= J2 a"''\^)PiPkPm+J2b'^'^{^)pmPg + W{^) = }C, + }C2 + ICo, {25) 

£,k,j,i=l m,q=l 
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This must satisfy the conditions 

2 

-2^a^"^(lnA), (26) 

s=l 
2 

-6 5^a^^^^(lnA)., z^j 

s=l 

2 2 

-^a^™(lnA),-3 5^a^'^'^(lnA)„ z^j (27) 

s=l s=l 
2 2 

6A^a^^'*V,-^6^^(lnA)„ (28) 

s=l s=l 
2 2 

2A J]a^**V,-^6^^(lnA)„ 

s=l s=l 
2 

Aj]6-v; = iy,. (29) 

Note that the a^'^-^* is a 4th order KiUing tensor. 

If is a 4th order symmetry operator, there exist functions a^^^^, W 
such that K has the unique self-adjoint form 

^ = E A^*^' (^''''^^'^'^) + (^''^^j) + = + ^2 + (30) 

where the functions b^^{xi, X2), Vr(xi, X2) contain the parameter dependence. 
Equating coefficients of the 5th derivative terms in the operator condition 
[K, H] = we obtain exactly the Killing tensor conditions 

The remaining conditions on K intertwine A, a^''^\ 6^*, W and V, and are 
complicated. Rather than solve them directly, we use the fact that the system 
with potential all must be solvable for all a, and require that the symmetry 
K[a) is analytic in a about a = 0. The following argument for the form of 
K is correct, up to addition of operators /(a)L2 or g{a) where L2 is a 2nd 
order self-adjoint a-independent symmetry operator. Modulo this remark, 
K[{a) must be at most a 1st order polynomial in a and -ft'o(a) must be at 
most quadratic. We can make the unique decomposition 



4af + af = 
3af" + 2af^ = 

21^ + bf = 
6f = 
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Substituting into [K{a), Ho + aU] = and equating the 3rd derivative terms 
that are hnear in a, we get exactly conditions ([2H1) Equating the coefficients 
of the 0th derivative terms that are quadratic in a we get exactly conditions 

Now we write K2{a) = A2 + OiB2,, K'q^q) = Aq + aBo + a'^Co. It follows 
that 

[K,, Ho] = 0, [K,, U] + [K2, Ho] = 0, [K2, U] + [Ko, Ho] = 0, (31) 
where 

K = K^ + K2 + Ko, Ki = K'^ + A2 + Ao, K2 = B2 + Bo, Ko = Co- 
Now define 

k = K,- K2U-\Ho + b)+Ko {U-\Ho + 6))' , 

where the operator order is important. A straightforward computation using 
identities, (13T|) yields [K, U~^{Ho + b)] = 0. Thus is a 4th order symmetry 
operator for the Hamiltonian U^^{Ho + b). 

Theorem 7 Let H(a) = Ho + aU, let K(a) be a 4th order self-adjoint 
symmetry of H{a), analytic at a = 0. Then there are 0th, 2nd and 4th 
order self-adjoint operators Kq, K2, ^4 such that K{a) = ^4 + aK2 + o^Kq 
and identities ( Tfil) are satisfied. The operator K = K4 — K2U~^{Ho + b) + 
Ko{U~^{Ho + &))^ is a 4th order symmetry for the system H = U~^{Ho + b). 

Corollary 8 K* = K so K is a 4th order formally self-adjoint symmetry of 
H. 

Example 5 This is an extension of Example\^to the quantum case, Let 
H = du+d22 + aixl + xl) + b-^^±^^+c^-^^^^ 

X2) X1X2 

There are two basic self-adjoint symmetry operators, one of 2nd order, 

22 / 4 I 4\ 

rf'-" ri^^ I rr>^ _l_ 1 

K2 = {x,d2 - X2d,f + + c^^ifl^ 

[xi — X — 2^y X\X2 



and one of 4th order. 



2 



K, = (dn - 022? + [2axl + 2b f^^^^ - 2c ^''\ l'" ]d-_ 



X2 /y>2\2 /7i2^2 
1 X2) X^X2 
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+ {2axi - ^)di + {2ax2 - ^)d2 + a^{xl - xlf + 



■^j^i. I v-""-^ ™3^^^ ' V-1 •"2J ' („2_„2\2 
Xl X2 \Xi X2) 



ii^! dl! + Sab ^^^^ + 2-^ + 6c(— + — ' 

12 \1 2/ tt'-^tju'^ jj-^ 



These operators generate a closed symmetry algebra. Let R — [K2, K^l.The 
relations are 

[K2, R] = 32H^K2 - S2{K4, K2} - 64(6 + 2c + 4)^4 + 64(6 - c + 2)//^ 

-128a(6 + l)Ki - 128a(6^ + 26c + 46 + 6c + 4), 
[Ki, R] = 32Kl-32H^K4+128aK2H^ -384:0^ Kl+128a{b+l)K4-64:a{b+4:C+6)H^ 

+256a2(2c - 6 + U)K2 + 128a^{b'^ + Ac^ + 206c + 186 + 8c - 8) 
There is also a Casimir constraint. Then the Stdckel transformed system 

H = T^, — T^ T^, — 2\ (^11 ~l~ ^22 + 

+ ^2) + ^^!2'^!2i + VJ'^ + ^ 

is also superintegrable with 4th and 2nd order generating self-adjoint symme- 
tries. 

5.4 Nth order operator symmetries 

A possible structure of the A^'th order operator case is now clear, though it is 
far from clear whether the structure inchidcs all cases. Suppose the system 
H{a) = Hq + V + aU = H + alf admits a truly A^th order symmetry operator 
K{a) analytic in a about a = 0, where N > 2 and K is self-adjoint for even 
N, skew-adjoint for odd N. Then we can write 



V2] 

^'N-2jy 



where each i^Ar-2j self-adjoint or skew-adjoint, depending on the parity of 
N. The symmetry condition is 

[N/2] 

[K + E ^^v-2,-(«), H + aU]^ 0, (32) 
3=1 
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where the K'j^_2j{a) are analytic in a. Suppose, modulo terms of the form 
f^'^\a)L]\f_2j where Ln-2j is an a-independent symmetry oi H + alJ for j > 0, 
we have 

= E ^!^-2,«\ J = 0, 1, ■ ■ ■ , [N/2] 

where the ^7v-2j independent of a. Setting KN_2j = X]L^=o~^''^^^^ ^!v-2j-2h 
we have K = YI^^q^ KN-2j and the symmetry condition (l32l) becomes 

[7V/2] 

a^irjv-2j, H + aU] = 0, (33) 

i=o 

or 

[KN-2jM + [KN^2j~2,H] = Q, j=0,l,--- ,[iV/2], (34) 

where we define -ft'iv_2j = for j > [N/2] and j < 0. 
Now define 

[N/2] 
h=0 

where 6 is a constant. From relations (1341) we have 

lN/2] 



[K,U-\H + b)] = J2i-^)''[KN-2h,U-\H + b)] {U-\H + b))' = (35) 

h=0 



[N/2] 

J2 {[KN-2h, U]{-1)^+\U-\H + b))^^' + [K^.2k, H]{-l)\U-\H + b)f) 

h=0 

[N/2] 

= $^(-1)'^ {[K^-2H-,2,U] + [K.2H,H]) {U-\H + b)f = 0. 

h=l 

Thus [k, H] = 0. 

Theorem 8 Let H{a) = Hq + V + aU = H + aU and N > 2. Let K{a) 
be a nonzero Nth order operator symmetry of H{a) analytic at a = 0, self- 
adjoint for even N and skew-adjoint for odd N. Suppose further there are 
operators Kj^_2j such that K(a) = Yl'jl/o^ -^JV-2jCt"' ond identities (E^j? i"^) 
are satisfied. Then the operator K = Y^\^Q\—l)'^KN^2h{U~^{H-[-b))^ is an 
Nth order symmetry for the system H = U^^{H + b). 
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Example 6 (The 9-1 anisotropic oscillator) This is a generalization of Ex- 
ample^ to a full Stdckel transform. Let H{0) = du + 822 + a{9xl + x^) and 
L be as in Example |^ with a replaced by a + a, and U = {9x1 + ^2) + ^- 
follows that the system 

H = ^g^2 /^2) ^ , (^11 + 922 + a{9x' + y') + b)) 
is superintegrable with one 2nd and one 3rd order symmetry. 

Note that Theorem [8] does not require that the quantum system go to a 
classical system, only that a scalable potential term can be split off. Thus it 
applies to "hybrid" quantum systems that have a classical part. 

Example 7 (The hybrid 9-1 anisotropic oscillator) Let H{0) = du + 822 + 
a{9xl + X2) — 2/x\. This is a superintegrable system with generating 2nd and 
3rd order symmetries, 

L = d22+axl, K = {xid2-X2di,d22} + {7:xl+ — ,di}-{3xi{axl-^^),d2}. 

6 X2 X2 

Note that this system does not have a classical limit. [Using a different nor- 
malization that makes clear the classical limit, Gravel writes this Hamiltonian 
as H{0) = -{hy2){dn + d22)+a{9xl + xl) + hyxlj Let U = {9xl + xl) + c. 
It follows that the system 

H = ,n 2 / , ( + ^22 + a{9xl + xl)-^ + b] 
[9x1 + 3;2) + c V ^2 / 

is superintegrable with one 2nd and one 3rd order symmetry. 



Example 8 (A translated hybrid 9-1 anisotropic oscillator) This is a slight 
modification of Example Let H{fS) = du + 822 + a{9xl + X2) + cxi — 
2/x\. This is a superintegrable system with generating 2nd and 3rd order 
symmetries, and no classical limit. Let U = Xi. It follows that the system 

H = {du + 822 + a{9xl + xl) + " + ^) 

is superintegrable with one 2nd and one 3rd order symmetry. This space is a 
Darboux space of type 1, 
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Example 9 Let H{0) = du + 822 + a/ x\ — 2/ x\. This is a superintegrable 
system with two linearly independent 2nd and three linearly independent 3rd 
order symmetries, J4I This system does not have a classical limit. [Using a 
different normalization that makes clear the classical limit, Gravel writes this 
Hamiltonian as H{0) = —{h^ /2){dii + 822) + a/ x\ + h^ / xl-] LetU = l/xj + c. 
It follows that the system 

H = (dii + d22 + ^ - ^ + b) 

1 + cxf \ xf X2 J 

is superintegrable with with two linearly independent 2nd and three linearly 
independent 3rd order symmetries. In the case c = this is a superintegrable 
system on a space of nonzero constant curvature. Indeed, for Xi,X2 real, it 
is the upper half space metric of non-Euclidean geometry. 

In the operator case where = in Theorem [8] there is always a corre- 
sponding classical system. Indeed, equations ( 133|) and ( IMl) clarify the close 
relationship between symmetries of quantum systems with potentials invari- 
ant under scaling and classical constants of the motion. If we set a = 1/^^, 
= in (IHHl) we can rewrite this expression as 



[iV/2] 



j=0 



Further if we write the differential terms in the operators Ki^_2j as 
KN-2j = ^ a*^'"*^"^^9ij ■ ■ ■dij^_^. + lower order terms, 

we can associate these operators with the phase space functions 

Then by equating coefficients of the highest order derivative terms in equa- 
tions (IMI) we obtain the Poisson bracket relations 

{/C^-2„ U} + {/C^_2,_2, ^0} = 0, J = 0, 1, ■ ■ ■ , [iV/2], (36) 

so that /C = Ei=o" _2j is an A^th order constant of the motion for the 
system with Hamiltonian Ti = Tio + U . 
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6 Conclusions and outlook 



We have found specializations of classical CCM that preserve the order of 
symmetries and determine symmetry algebra homomorphisms, and for 2D 
manifolds we have extended them to the quantum case. Generally speaking, 
these transforms apply to systems with a nonconstant potential that admits 
scaling in at least one parameter. They do not apply to quantum systems 
with no classical counterpart in which the potential is fixed. This tool makes 
it clear that superintegrable systems occur for a wide variety of manifolds, 
not just on constant curvature spaces. For 2nd order superintegrable systems 
the Stackel transform has been used effectively in 2D to show that all such 
systems are transforms of systems on constant curvature spaces, and this 
has lead to an elegant classification of all such systems. It is our aim to 
develop CCM to investigate the possibility of a similar classification for 3rd 
and higher order superintegrable systems. 

For simplicity, we have restricted our quantum constructions to 2D man- 
ifolds though some partial results hold in n dimensions. There appears to 
be no insurmountable barrier to extending these results to 3D and higher 
conformally fiat manifolds, hut the details have not yet been worked out. 
Clearly gauge transformations are required and the gauge will be a function 
of the scalar curvature of the manifold. 
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